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O . ABSTRACT 

' Aims. We describe the mid- (MIR) and far- (FIR) infrared properties of a large (~1000) sample of the most isolated galaxies in the local 
II Universe. This sample is intended as a "nurture-free" zero point against which more environmentally influenced samples can be compared. 
Q Methods. We reprocess IRAS MIR/FIR survey data using the ADDSCAN/SCANPI utility for 1030 out of 1050 galaxies from the Catalogue 
^ of Isolated Galaxies (CIG) as part of the AMIGA project. We focus on diagnostics (FIR luminosity LpiR, R = \og{Lpifi/ Lg) and IRAS colours) 
C/2 thought to be sensitive to effects of environment or interaction. 

Results. The distribution of log(LFiR) is sharply peaked from 9.0-10.5 with very few (<2%) galaxies above 10.5. Review of available optical 
^ , images of the most FIR luminous galaxies finds the majority to be, likely, interacting systems missed in our earlier morphological reevaluation. 
T ""j ■ The optically normalised luminosity diagnostic R = log(Lj:jR/ Lg) shows a distribution sharply peaked between 0.0 and -1.0. These results 
I^nJ ' were compared to the magnitude limited CfA sample that was selected without environmental discrimination. This modestly (e.g. compared to 
cluster, binary galaxy and compact group samples) environmentally affected sample shows significantly higher mean log(LpiR) and R, whereas 
the mean log(LB) is the same. Our sample shows a strong LpiR vs. Lg correlation, with a slope steeper than one (LpiR ^b"*')- Interacting 
galaxies were found above this correlation, showing an enhancement in Lfir. With respect to the IRAS colours, we found higher Fgo/Fjoo value 
for ellipticals and late-type galaxies than for spirals, indicating a higher dust temperature. The mean value of -FjoZ-Fioo was found to be lower 
than for interacting samples from the literature. 

Conclusions. The results indicate that the FIR emission is a variable enhanced by interaction, and that our sample probably shows the lowest 
possible mean value. This attests to the utility of our sample for defining a nurture-free zero point. 

Key words, galaxies: evolution - galaxies: interactions - galaxies: luminosity function, mass function - galaxies: ISM - surveys - infrared: 
galaxies 



1. Introduction 

Although it is widely accepted that galaxy interactions can 
stimulate secular evolutionary effects in galaxies (e.g. en- 
hanced star formation, morphological peculiarities including 
transitions to earlier type, active nuclei) (e.g. Sulentic 19761 
Hernquist'T989', Xu & Sulentic 1991 1 there are still many open 
questions. Studies aimed at quantifying the level of interaction 

* Full Tables|2|and|3|are available in electronic form at the CDS via 
anonymousftptocdsarc.u-strasbg.fr (138.79.128.5) or via 
http : // cdsweb . u-strasbg . fr/ cgi-bin/q cat?3/A+A/vvv/ppp, 
and from http : //www . iaa . es/AHIGA. html. 



enhancement have even produced contradictory results; e.g. 
some studies of interacting pairs find a clear star formation en- 
hancement (Bushouse et al. 119871 Bushouse[1988 ) while others 
find only a marginal increase (Bergvall et al. lZ003t . Much of 
this uncertainty reflects the lack of a statistically useful base- 
line. What is the amplitude and dispersion in a given galaxy 
property that can be ascribed to "nature"? A definition of "iso- 
lated galaxy" is needed before one can properly assess the his- 
tory and properties of non-isolated ones. This has motivated us 
to build a well-defined and statistically significant sample of 
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isolated galaxies to serve as a control sample for the study of 
galaxies in denser environments. 

The AMIGA project (Analysis of the interstellar Medium 
of Isolated GAlaxies) involves the identification and study of 
a statistically significant sample of the most isolated galax- 
ies in the local Universe. Our goal is to quantify the prop- 
erties of different phases of the interstellar medium in these 
galaxies which are likely to be least affected by their external 
environment. We adopted the Catalogue of Isolated Galaxies 
(CIG: Karachentseva 1973, Karachentseva et al. 1986l, includ- 
ing 1051 objects, as a base sample. All CIG galaxies are part 
of the Catalogue of Galaxies and Clusters of Galaxies provid- 
ing reasonably uniform apparent magnitude measures (CGCG: 
Zwicky et al. 11961-19681 with m^g < 15.7 and 6 > -3deg. 
Redshifts are now virtually complete for this sample with only 
one of the compiled objects recognised as a Galactic source 
(CIG 781 = Pal 15; Nilson [T973t reducing the working sample 
to « = 1050 objects. AMIGA is compiling data that will char- 
acterise all phases of the ISM: blue magnitude, mid- and far- 
infrared. Ha, and radio continuum fluxes, as well as the emis- 
sion of the atomic gas (HI) and of carbon monoxide (CO), as 
a tracer of the molecular gas. The data are being released and 
periodically updated at http : //www . iaa . es/AMIGA . html 

Previous AMIGA papers evaluated, refined and improved 
the sample in different ways including: 1) revised positions 
(Leon & Verdes-Montenegro 2003), 2) sample redefinition, 
magnitude correction and full-sample analysis of the Optical 
Luminosity Function (OLF) (Verdes-Montenegro et al. '2005' 
Paper I) and 3) morphological revision and type-specific OLF 
analysis (Sulentic et al. 120061 Paper II). In the present paper 
we analyse basic mid- (MIR) and far-infrared (FIR) properties 
using data from the IRAS survey (Sects. |2] and O. In Sect. |3 
of the paper we present the FIR luminosity function followed 
by consideration of various MIR and FIR diagnostics that have 
been used in the past to quantify the effects of environment. In 
Sect.|5lwe discuss the results and compare them to other stud- 
ies. Future papers will present a quantification of the isolation 
condition and the analysis of the radio continuum. Ha, CO and 
HI data. 

2. ADDSCAN/SCANPI analysis of the IRAS data 

We present co-added ADDSCAN/SCANPI derived fluxes or 
upper limits for 1030 AMIGA galaxies. The remaining 20 
galaxies in our sample were not covered by the IRAS survey. 
Previous studies involving CIG galaxies worked with smaller 
subsamples and, in most cases, used IRAS data from the Point 
Source (PSC) and Faint Source Catalogues (FSC). A subsam- 
ple of 476 CIG galaxies with redshifts and PSC fluxes were 
used as a control sample for a study of FIR emission from 
isolated pairs (Xu & Sulentic [T99T1 hereafter XS91). Verdes- 
Montenegro et al. ( 1998 1 constructed a reference sample of 68 
CIG galaxies with redshift and blue luminosity distributions 
matching their target set of Hickson (il98 2i compact groups. 
Hernandez-Toledo et al. ( fT999t obtained SCANPI data for 465 
CIG galaxies (those with available redshift data) in order to use 
them as a reference in a study of galaxy pairs. FIR data for the 
CIG galaxies were however not published in that paper 



Table 1. Detection rates and point vs. extended source numbers 
for the CIG IRAS sample (n = 1030). 



A 


Threshold 


Detections 


Detection rate 


Extended 


Point 


12 


30- 


180 


17% 


40 


141 




50- 


94 


9% 


37 


57 


25 


30- 


245 


24% 


57 


188 




50- 


158 


15% 


53 


105 


60 


30- 


729 


71% 


84 


645 




50- 


591 


57% 


82 


509 


100 


3cr 
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IRAS PSC and FSC data exist for only about half of the 
galaxies in our sample motivating us to ADDSCAN/SCANPI 
reprocess our entire sample. We used the revised posi- 
tions from Leon & Verdes-Montenegro (12003> which have 
a precision of 0'.'5. ADDSCAN/SCANPI, a utility provided 
by the Infrared Processing and Analysis Center (IPAC) 
(http://scanpi.ipac. caltech.edu: 9000/1, is a one- 
dimensional tool that coadds calibrated IRAS survey data. It 
makes use of all scans that passed over a specific position and 
produces a scan spectrum along the average scan direction. It is 
3-5 times more sensitive than IRAS PSC since it combines all 
survey data (Helou et al. 1988 1 and is therefore more suitable 
for detection of the total flux from slightly extended objects. 
Our sample was well suited for ADDSCAN/SCANPI process- 
ing because: i) confusion is minimised since our sample galax- 
ies were selected with an isolation criterion and ii) the galax- 
ies are small enough to permit derivation of reliable fluxes. 
An analysis of the IRAS Bright Galaxy Sample (BGS) with 
ADDSCAN/SCANPI (Sanders et al. 12003^ found that missed 
flux became important only for optical sizes larger than 25'. 
About 97% of the galaxies in our sample are smaller than 4'. 

ADDSCAN/SCANPI derives four different flux estimators: 
a) Peak: maximum flux density within the signal range speci- 
fied, h) friu(z): total flux density estimated from integration of 
the averaged scan between the zero crossings, c) Tempi: flux 
density estimated from the best-fitting point source template 
and d) fnuit): total flux density estimated from integration of 
the averaged scan between fixed points defining an integration 
range. We adopted the default SCANPI ranges (correspond- 
ing to the nominal IRAS detector size) [-2', -1-2'], [-2', +2'], 
[-2:5, -h2:5] and [-4', H-4'] at 12, 25, 60 and 100 //m, respec- 
tively. We used the median as the most robust combination of 
scans and followed IPAC indications in order to choose the best 
flux density from among the estimators for each galaxy. We 
first flagged as detected those galaxies with a S /N > 3. We 
visually confirmed all cases and found some errors produced 
by bright stars in the field or baseline corruption from noise or 
cirrus. 

Table^summarises the number of detected sources at each 
IRAS band. For completeness, we also include the correspond- 
ing numbers for a detection threshold of 5 times the rms noise 
level, which is the limit used in the data analysis carried out 
in this paper (see Sect. 14. 1> . Figure^shows the detection rate 
(at a 3o- detection threshold) at the four IRAS wavelengths as 



Lisenfeld et al.: The AMIGA sample of isolated galaxies III 



3 




-5 5 10 

Hubble type 



Fig. 1. Detection rate at 12, 25, 60 and 100 /im as a function of 
Hubble type. 

a function of Hubble type. The MIR-FIR detection rates show 
a minimum for early-type galaxies gradually increasing from 
10-20% to 20-80% for late-type spirals. We see a decline to 
20-60% beyond type Sd (T = 7) probably reflecting an in- 
creasing dwarf galaxy population with low dust masses. 

Figure|2lplots the "miss" parameter which is the off'set in ar- 
cmin between the galaxy position and the position of the signal 
peak along the average scan direction. This parameter is used 
as the primary measure of source identification. The majority 
of sources cluster around zero offset with the largest deviations 
occurring at 12 and 25 //m because: 1) the resolution is higher 
and 2) the S /N is usually lower than at longer wavelengths. 
The standard deviations of the "miss" parameter are 18", 24", 
14" and 28", respectively, for 12, 25, 60 and 100 fim. This is a 
factor of 5 to i of the nominal FWHM of the IRAS detectors 
(0:77, 0:78, 1:44, and 2:94 at 12, 25, 60 and 100 jum respec- 
tively). This scatter is reasonable when one allows for the fact 
that most of these galaxies are not very infrared (IR) bright so 
that determination of the source centroid depends somewhat on 
the S IN of the measurement. 

In the next step we used two different tests to decide 
whether a detected source was extended or pointlike with re- 
spect to the IRAS beam. In Test 1 we considered as extended 
those galaxies where the signal width was greater than the ex- 
pected width for a point source. We used both w25 and \v50 
(width of the signal in arcminutes at 25% and 50% peak) for 
this comparison. We compared our measures to the widths of 
point sources (Sanders et al. I2003t where w25psf wSOpsf 
were 1.'40, 1.'38, 2.'06, 432 and 1.'04, i;00, i;52, 3:22 at 
12, 25, 60, 100 jum, respectively. In Test 2 we classified as 
extended those sources where the integrated flux, fnu(z) was 
substantially larger than the peak flux adopting the condition 
fau{z) - Peak > 5cr as a threshold criterion for an extended 
source. The percentages of galaxies showing conflicting clas- 



sifications in the two tests were 9, 17, 23 and 18% at 12, 25, 
60 and 100 fim, respectively. We revised these cases interac- 
tively and found the differences were most often due to baseline 
corruption by noise and/or cirrus. Tabled lists the number of 
sources classified as point or extended for each IRAS band. The 
5cr cutoff reduces, compared to the 3cr cutoff, mainly the num- 
ber of detected point sources and leaves the number of detected 
extended sources almost unchanged. The reason is that sources 
classified as extended have in most cases S /N > 5. Once the 
size of a source was decided we could choose a flux estimator 
following guidelines given by IPAC. For point sources three 
cases were considered: a) if Tempi > 20 Jy we used Peak, b) 
if Tempi < 1 Jy (< 2 Jy at 100 pm) we used Tempi and c) if 
1 Jy < Tempi < 20 Jy (2 Jy < Tempi < 20 Jy at 100 yum) we 
used Tempi if Peak and fnu{t) agreed within 3cr. Otherwise we 
visually determined the best estimator of the total flux density. 
In the case of extended sources we used/nM(z) when S /N > 10 
and/«M(f) for the rest. 

3. The data 

Table|2llists A 12, 25, 60 and 100 pm derived fluxes obtained 
using procedures explained in Sect. |2l We also tabulate some 
related parameters, as detailed in the notes to the table. 

3.1. Comparison to other IRAS catalogues 

We compared IRAS fluxes obtained with SCANPI to data 
available from the IPAC archives and in the literature. We re- 
trieved data from the Faint Source Catalogue (FSC) and the 
Point Source Catalogue (PSC) from the IRAS database through 
the GATOR service (http://irsa.ipac.caltech.edu/i. 
We found 509 CIG galaxies in the FSC and additional data 
for 15 galaxies in the PSC. The average error-weighted ra- 
tios F(SCANPI)/F(FSC+PSC) for galaxies detected both by 
SCANPI and in the FSCh-PSC are 1.24 + 0.50 (n = 114), 
1.16 + 0.33 (n = 153), 1.09 + 0.15 (n = 501), 1.05 + 0.13 
(n = 407) for 12, 25, 60 and 100 pm. The average ratio is 
slightly larger than one and decreases with increasing wave- 
length. This indicates that the flux derivation with SCANPI is 
able to pick up more flux for extended objects than FSC/PSC 
especially at short wavelengths were the IRAS beam is smaller. 
There is a large number of galaxies with FSC/PSC tabulated 
upper-limits (n = 55, 70, 9, and 81 for 12, 25, 60 and 100 pm) 
that were replaced by SCANPI detections indicating that the 
detection rate has been improved by the reprocessing. Other 
galaxies were listed as FSC/PSC detections while SCANPI de- 
rived only upper limits (n = 29, 21, 5, and 3 for 12, 25, 60 and 
100 pm). We checked those cases individually and found that 
all were weak sources where either: 1) the automated SCANPI 
procedure did not confirm a detection or 2) we decided, after 
visual inspection, that 5 /A^ < 3. 

Figure |3] shows the SCANPI-to-FSC/PSC flux ratio as a 
function of optical diameter for each IRAS band. We ex- 
clude detections where SCANPI and FSC/PSC fluxes agree 
within the uncertainties (which is the case for 82, 107, 397 
and 368 galaxies at 12, 25, 60 and 100 pm, respectively). 
The 12 and 25 pm plots, and to a lesser extent 60 pm, show 
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Fig. 2. Histogram of the "miss" parameter (offset in arcmin between the galaxy position and the position of the signal peak) for 
each IRAS band. 



Table 2. FIR flux densities ^ 
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0.12 


4 


y 


5 


<0.13 
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0.12 
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The Table format is: Column I: CIG number. Column 2: Flux density at 12 yum, calculated as explained in Sect.|2] Upper limits are preceded 
by a "<" sign. A 3(T value has been adopted for the upper limits, except for CIG 397 where the 12 /im scan presents confusion with a close star 
and 20% of the peak of the emission has been adopted as an upper limit. Column 3: rms noise of the data at 12 fim. Column 4: Method used 
to derive the flux densities given in column (2). Codes 1 to 4 correspond to the following flux estimators: 1 = Tempi, 2 = Peak, 3 =fnu(t), 4 = 
fnu(z). Code 5 corresponds to upper limits obtained as 3cr. Code 6 is reserved for some particular cases: CIG 397 (see above) and nine galaxies 
included in the catalogue of large optical galaxies of Rice et al. 1 1988 1 (CIG 105, 197, 324, 347, 461, 469, 523, 559 and 610) where we have 
used the values given in their catalogue (see also Sect. l3.lt . Column 5: Detected galaxies are flagged with "y" if they have been classified as 
extended, and with "n" when classified as point sources. Column 6-9: The same as column 2-5 for 25 yum. Column 10-13: The same as column 
2~5 for 60 ytim. Column 14-17: The same as column 2-5 for 100 jim. 

' The full table is available in electronic form at the CDS and from fhttp : //www . iaa . es/AMIGA . html | 



that F(SCANPI)/F(FSC+PSC) increases with optical diameter 
above about l'-3'. This supports our inference that FSC/PSC 
flux densities are often underestimated for large galaxies be- 
cause part of the flux falls outside the IRAS beam and that 



SCANPI is able to provide a more realistic flux estimate for 
these sources. There are only three galaxies (CIG 546, 616 
and 721) with F(SCANPI)/f (FSC+PSC) significantly below 
one. In the case of CIG 546 (F6o(SCANPI)/F6o(FSC+PSC) 
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Fig. 3. SCANPI/FSC+PSC flux ratio as a function of optical diameter for the 4 IRAS wavelengths. Detected galaxies showing 
consistent fluxes between SCANPI and FSC+PSC are not included in the figure for clarity. 



= 0.66) and CIG 721 (F6o(SCANPI)/F6o(FSC+PSC) = 0.72) 
the origin of the difference is unclear since we have no 
reason to doubt the reliability of our SCANPI estimates. 
CIG 616 (Fi2(SCANPI)/Fi2(FSC+PSC) = 0.56) was detected 
by SCANPI just above a 3cr level so that the flux density has a 
larger uncertainty. 

We compared 18 galaxies in common with the Bright IRAS 
Galaxy sample (Feo > 5.24 Jy) where flux densities were also 
derived using SCANPI (Sanders et al. 2003 1. Agreement is bet- 
ter than 6% for all sources at 12, 60 and 100 /urn. At 25 yum 
there are three sources (CIG 442, 549 and 1004) where our 
adopted values exceed those derived by Sanders et al. (I2003> at 
the 10-20% level. We think that the discrepancy arises because 
some of the flux in these sources extends beyond the integra- 



tion range used in deriving /««(?) and will therefore be better 
estimated using our/«M(z) values. 

Following IPAC recommendations we compared the re- 
sults derived with SCANPI to those derived from 2D Full 
Resolution Coadded (FRESCO) images for sources with op- 
tical diameter larger than 2^5 (107 objects). Since FRESCO 
images do not have large-scale background removed (they 
are not point-source filtered), they provide additional infor- 
mation about galaxy environments including possible confu- 
sion due to nearby stars or Galactic cirrus. We extracted in- 
dividual source fluxes from FRESCO images using Sextractor 
(Bertin & Arnouts 1996). We extracted fluxes for CIG galax- 
ies using both 3 and 5cr thresholds above the local back- 
ground level in order to estimate the effects of background 
and particularly Galactic cirrus. Calibrated FRESCO fluxes 
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for the 4 IRAS bands were compared to the SCANPI fluxes 
and we found ratios F(SCANPI)/F(FRESCO) close to unity 
(1.04+0.42, 0.98+0.40, 0.89+0.33 and 0.97+0.44 at 12, 25, 60 
and 100 fim) respectively. Scatter was high but we did not find a 
trend with optical diameter that might point towards flux being 
missed using either procedure. More likely contamination from 
the local foreground (Galactic emission) is dominating the flux 
determination. 

Finally, we searched for CIG galaxies included in the cata- 
logue of large optical galaxies (Rice et al. J.988 1 and found nine 
objects (CIG 105, 197, 324, 347, 461, 469, 523, 559 and 610). 
In most cases we find reasonable (within 2cr) agreement be- 
tween flux estimates. There are some significant discrepancies 
for the largest objects. The most severe discrepancy involves 
the galaxy with largest apparent optical diameter CIG 610 (= 
MlOl, 28^8 X 26:9) where SCANPI fluxes are only 10-35% of 
the fluxes given in Rice et al. J1988> . Disagreements of up to a 
factor of 2.5 are present for CIG 197 (= NGC 2403, 21:9x123) 
and CIG 523 (= NGC 4236, 21.9' x 7.2'). The single dis- 
crepancy involving a smaller galaxy, CIG 105 (= NGC 925, 
10:5 X 5:9) finds a SCANPI flux at 12 ;um that is almost a fac- 
tor two lower than the flux in Rice et al. ( 1988 1. We adopted 
the flux densities of Rice et al. J1988> for all nine galaxies. 

4. Data analysis 

4.1. Sample definition 

In the following sections we analyse the FIR emission proper- 
ties of the CIG galaxy sample. In order to do this in a statisti- 
cally meaningful way we focus on the optically complete sam- 
ple described in Paper I. This sample involves galaxies with 
corrected Zwicky magnitudes in the range 1 1 .0-15.0 for which 
we found < V/V„i >= 0.40, indicating 80% completeness. 

We include in the present work some changes/upgrades 
with respect to the previously (in Paper I) defined sample: 

1. We include 20 galaxies for which redshift informa- 
tion has become available since the publication of 
Paper I (the updated redshift list is available from 
http : //www . iaa . es/AMIGA . html). 

2. Morphological revision of the sample, described in Paper II 
identified 32 galaxies that are probably not isolated in 
the sense that they might involve isolated interacting pairs 
and/or multiplets. These galaxies are excluded from the 
most isolated sample and represent part of the AMIGA re- 
finement. However they provide us with a useful internal 
comparison sample to test the effects of interaction con- 
tamination. 

3. We recomputed corrections to the Zwicky magnitudes fol- 
lowing Paper I but using the revised morphologies from 
Paper II. This change in individual magnitudes will there- 
fore slightly change the sample involving galaxies in the 
range 1 1 .0-15.0 mag. The present sample shows a value of 
< V/V,„ >= 0.43 indicating a slightly improved level of 
completeness compared to Paper I. 

4. We exclude two nearby dwarf ellipticals (CIG 663 = Ursa 
Minor and CIG 802 = Draco) for which we have only 



Table 3. FIR and blue luminosities^ 



(1) 


(2) 


(3) 


(4) 


CIG 


Distance 


log(LFIR) 


log(Lfi) 




(Mpc) 




(Lq) 


1 


92.2 


10.23 


10.44 


2 


88.7 


< 9.72 


9.76 


3 
4 


26.1 


9.91 


10.17 


5 


100.2 


9.75 


10.07 



The entries are: Column (1): CIG number. Column (2): Distance in 
Mpc from Paper I. Column (3): FIR luminosity, calculated as de- 
scribed in Sect. 14.21 Upper limits are indicated with < in front of the 
value. Galaxies with distances, but without FIR data points (in total: 
20 objects) lie in the area not covered by IRAS. Column (4): Blue lu- 
minosity, calculated as described in Sect. l4.4l 

' The full table is available in electronic form at the CDS and from 
http : //www . iaa . es/AMIGA . html 



IRAS upper limits and very low inferred luminosity limits 
(log(LnR/io) < 3.25). 

We are left with a sample of 719 galaxies with IRAS 
data, and redshift data is available for 701 galaxies of them. 
Hereafter we will refer to this sample as the AMIGA (FIR) 
sample. We decided to increase the detection threshold to 5 
cr in order to make sure that we only consider reliable detec- 
tions. Thus, with respect to Table |2 we now consider only 
those fluxes as detections where the S /N ratio is above 5, and 
we use an upper limit of 5 times the rms noise for values be- 
low. (We chose to leave the 3cr detection limit in Table |2] in 
order to provide the complete data set.) With this higher thresh- 
old, 511 galaxies have a detection at least at one wavelength. 
This sample can be cut in many diff'erent ways. Right now 
we make no restriction in recession velocity. This allows us 
to sample the widest possible luminosity range. Sources with 
V < 1500 kms ' provide an insight into the IR emission from 
local dwarf galaxies that are not included in the rest/bulk of 
the sample. The drawback about including these galaxies in the 
sample involves the difficulty in reliably assessing their isola- 
tion properties. 

4.2. FIR luminosity 

FIR luminosity, Lfir, is computed from IRAS measurements 
as log(LFiR/i0) = logiFIR) + 21og(D) + 19.495, where D is 
distance in Mpc and = 1.26xlO"''*(2.58F6o + Fioo) Wm^^ 
(Helou et al. [T988t the flux in flie FIR range, with the IRAS 
fluxes at 60 and 100 /im, F(,o and Fioo. Lfir and the distances 
adopted are listed in Table |3] 

Figure 0] shows the distribution of FIR luminosity for the 
optically complete AMIGA sample and in Table 0] the mean 
and median values are given. We include in Fig. |4] individual 
histograms for: 1) galaxies detected at both 60 and 100 jjm, 
2) those not detected at one or both wavelengths and 3) the 
distribution calculated using survival analysis that takes upper 
limits into account. We use the ASURV package for the lat- 
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Table 4. Mean and median values of Lfir. I^b^ and R - ^og(L^YB.I Lg). 



1 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Sample 


n 


< log(Ls) > 




< l0g(Z.FiR) > 


n 




<R> 






med(log(LB)) 




med(log(LFiR)) 






med{R) 


Total 


701 


9.97±0.02 


312 


9.15+0.06 


719 


327 


-0.56+0.03 






10.06 




9.74 






-0.29 


S/Im {T = 1-10) 


616 


9.98+0.02 


248 


9.26+0.05 


634 


263 


-0.49+0.02 






10.07 




9.76 






-0.30 


E (T = -5) 


27 


9.95+0.06 


21 


8.83+0.16 


27 


21 


-1.01+0.12 






10.01 




9.77 






-0.14 


SO {T = -2) 


36 


9.82±0.08 


27 


8.58+0.15 


36 


27 


-0.95+0.09 






9.95 




9.65 






-0.23 


SOa (T = 0) 


14 


9.88+0.07 


8 


9.33+0.16 


14 


8 


-0.57+0.15 






9.88 




9.73 






-0.03 


Sa(r = 1) 


10 


9.92+0.18 


3 


9.27+0.20 


10 


3 


-0.60+0.10 






10.00 




9.64 






-0.49 


Sab (T = 2) 


39 


10.05+0.05 


15 


9.39+0.10 


39 


15 


-0.62+0.08 






10.00 




9.61 






-0.29 


Sb (T = 3) 


115 


10.06+0.04 


42 


9.54+0.08 


118 


45 


-0.38+0.04 






10.10 




9.90 






-0.14 


Sbc (T = 4) 


155 


10.10+0.03 


70 


9.37+0.09 


160 


73 


-0.50+0.03 






10.15 




9.82 






-0.32 


Sc (T = 5) 


182 


10.12+0.03 


69 


9.62+0.05 


188 


74 


-0.38+0.02 






10.20 




9.91 






-0.28 


Scd (T = 6) 


47 


9.65+0.08 


19 


8.89+0.12 


47 


19 


-0.54+0.06 






9.77 




9.28 






-0.35 


Sd (T = 7) 


34 


9.58+0.09 


15 


8.73+0.15 


38 


19 


-0.64+0.06 






9.51 




8.95 






-0.52 


Sdm(r = 8) 


10 


9.38+0.21 


7 


8.40+0.19 


10 


7 


-0.55+0.05 






9.05 




8.49 






-0.50 


Sm (T = 9) 


9 


9.11+0.34 


5 


7.98+0.27 


9 


5 


-0.58+0.05 






9.07 




8.56 






-0.61 


Im(r = 10) 


15 


9.01+0.21 


3 


8.30+0.34 


15 


3 


-0.58+0.13 






9.04 




8.70 






-0.47 


Interacting 


14 


9.99+0.11 


2 


9.87+0.20 


14 


2 


-0.06+0.08 






9.98 




10.02 






-0.11 



The entries are: Column 1: Subsample considered. All subsamples are selected from the optically complete, magnitude limited subsample. The 
interacting subsample consists of the galaxies from the CIG excluded in Paper II (see Sect. l4.lt . Column 2: Total number of galaxies with 
velocity and IRAS data in the subsample. Column 3: First row: Mean value of Lg. Second row: Median value of Lg. Column 4: Number of 
upper limits in FIR (at 60 or 100 fim). Column 5: First row: Mean value of Lfir, using the Kaplan-Maier estimator from ASURV. Second row: 
Median value of LpiR, only for detections. Column 6: Total number of galaxies with IRAS data in the subsample. Column 7: Number of upper 
limits in FIR (at 60 or 100 ^m). Column 8: First row: Mean value of R = log(Lpis,/LB), using the Kaplan-Maier estimator from ASURV. Second 
row: Median value of R, only for detections. 



ter calculations throughout this paper'. The distribution peaks 
in the bin log(LFiR/i0) = 9.5-9.75 with the ASURV esti- 
mated mean log(LFiR/i'o) - 9.15 (see Table 0}. Practically 
all galaxies have FIR luminosities between log(LFiR/io) = 7.5 
and log(LFiR/Lo) - 11-25. Only one object, the faint irregu- 
lar member of the Local Group CIG 388 (= Sextans B), shows 
log(LFiR/Lo) = 6.01. It is remarkable that the bulk of the FIR 
luminosities (98%) lies below log(LFiR/L0) = 10.5. 

In Fig.|5lwe show the distribution of the FIR luminosities 
as a function of Hubble type, and in Table |3 we list the mean 
(taking upper limits into account with ASURV) and median 



' Astronomy Survival Analysis (ASURV) Rev. 1.1 is a gener- 
alised statistical package that implements the methods presented by 
Feigelson & Nelson 1,19851 and Isobe et al. U986i) . and is described in 
detail in Isobe & Feigelson <1990t and La Valley et al. <1992L 



(for detections only) values. The mean values show a strong 
increase after the early-types (E-SO) beginning at T = (SOa) 
and increasing through T = 5 (Sc), followed by a decline to 
a near constant mean for the latest types {T = 8-10). ASURV 
means for early-types are most strongly driven by upper limits 
with most detected E-SO showing Lfir values above the com- 
puted means, similar to those for late-type spirals. This marks 
the detections as unusual indicating that these may not be typ- 
ical (or even) E-SO galaxies (see discussion in Paper II). As 
we proceed from left to right in the plot the effect of upper 
limits gradually decreases and mean and median values con- 
verge. Our previously identified (Paper II) dominant (~65%) 
isolated late-type (T = 3-5) spiral population shows FIR lu- 
minosities strongly concentrated (due to the minimisation of 
nurture-driven dispersion) in the range 9.4-10.5. We also ob- 
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50 



loc(Lp,„) [LJ 

Fig. 4. The FIR luminosity distribution of the optically selected 
sample described in Sect. 14. II The full line shows the distribu- 
tion calculated with ASURV (n - 701), the shaded area shows 
all galaxies detected at both 60 and 100 yum (n = 478), and the 
dashed line gives the non-detections. 




Hubble type 

Fig. 5. Distribution of FIR luminosity as a function of Hubble 
type. Only detections are shown. The open triangles give the 
mean value for each Hubble type, calculated with ASURV and 
taking the upper limits into account. The open squares are the 
median values for the detections only. 



serve a small but significant population of spiral types T - 2- 
7 with very low FIR luminosities. We see an apparent strong 
drop in mean FIR luminosity (~0.7 in log(LpiR)) later than type 
r = 5. If real, there are three candidate explanations: 1) de- 
creasing dust mass simply following decreasing galaxy mass 



for Scd-Sd, 2) decreasing dust content in Scd-Sd or 3) less ef- 
ficient star formation in Scd-Sd (always relative to Sb-Sc). The 
latest types show minimal upper limits since they are very local. 
This mostly dwarf galaxy population falls out of our magnitude 
limited sample beyond a few xlO^ kms"' recession velocity. 

4.3. FIR luminosity function 

Since the AMIGA sample is optically selected we derive the 
FIR luminosity function (FIRLF) from the optical luminosity 
function and the fractional bivariate function between FIR lu- 
minosity and optical luminosity (see Paper I). The differential 
FIRLF, which gives the number density of galaxies per unit 
volume and per unit log(LpiR) interval is derived from the fol- 
lowing formula: 



^(L) = 2.5AM y @(L\Mj)^(Mj), 



(1) 



where L = log(LFiR), and is the differential FIRLF. The bi- 
variate (optical, FIR) luminosity function @{L\Mj) is defined 
as 



&(L\Mj) 



P{L,Mj) 
AL ' 



(2) 



where AL = 0.25 and P{L,Mj) is the conditional probabil- 
ity for a source with absolute magnitude M (Mj + 0.5AM > 
M > Mj - 0.5AM) to have the logarithm of its FIR luminos- 
ity, log(LFiR), within the interval [L - 0.5AL,L + 0.5AL]. The 
Kaplan-Meier estimator (Schmitt 1985; Feigelson & Nelson 
[1985^, which also exploits the information content of upper 
hmits, has been used in computing the bivariate luminosity 
function and the associated errors. <i> is the differential OLF 
per unit volume and per unit magnitude interval, AM is the bin 
width of the OLF in magnitude units. The factor 2.5 arises be- 
cause a unit magnitude interval corresponds to only 0.4 unit 
of L. The summation is over all bins of the OLF. The errors of 
*F(L) are the quadratic sum of the uncertainties for the OLF and 
bivariate luminosity function. 

Figure |6l shows the resultant FIRLF and Table |5l lists the 
corresponding values. Also shown is the FIRLF from XS91 for 
a smaller subsample of the CIG. We see that the general shape 
has not changed substantially with the use of a larger and more 
complete sample. It is our contention that it represents the best 
"natural" or "nurture-free" FIRLF yet derived. A strong decline 
in the FIRLF above log(LFiR/L0) ~ 10.5 is clearly visible. In 
the last few bins there are only very few objects (see Table|5}, 
making the value of the FIRLF uncertain. 

We have fitted the FIRLF with a Schechter function: 



^(L).^o(^fexp(-A). 



(3) 



The best-fit parameters are - (7.4 + 1.4) x lO"** Mpc"^ 
(6log(LnR)r\L* = (1.9+O.2)xlOi"L0 and 0- = -0.46+0.05. 
We have also plotted in Fig.|5]the fit to the IRAS Bright galaxy 
Sample FIRLF (Sanders et al. 2003). They found, in agreement 
with other FIR selected samples, that a double-power law pro- 
vides the best fit to the data. The difference from a Schechter 
fit typically starts to be noticeable above 1O'^L0. Sulentic & 
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Fig. 6. Bivariate FIR luminosity function of our sample com- 
pared to the CIG sample used in XS91. The dotted line is the 
fit with a double power-law derived in Sanders et al. ('2003 1 for 
the Bright IRAS galaxy sample. The full line is a Schechter fit 
to our data. 

Table 5. FIR luminosity function 



log(LFiR) (in Lq) y(Mpc-3log(L)-') n 



7.50 


4.15 X 10- 


-3 


± 2.82 X 10- 


-3 


2.5 


7.75 


4.15 X 10 


-3 


+ 2.82 X 10- 


-3 


2.5 


8.00 


1.09 X 10 


-2 


+ 3.05 X 10- 


-3 


19.6 


8.25 


6.09 X 10- 


-3 


+ 2.17X 10- 


-3 


28.7 


8.50 


4.85 X 10- 


-3 


+ 1.32X 10- 


-3 


23.0 


8.75 


3.73 X 10- 


-3 


+ 1.07X 10- 


-3 


34.9 


9.00 


3.01 X 10- 


-3 


± 5.63 X 10- 


-4 


101.2 


9.25 


1.02 X 10- 


-3 


+ 2.35 X 10- 


-4 


38.1 


9.50 


1.81 X 10- 


-3 


+ 2.56 X 10- 


-4 


123.9 


9.75 


1.25 X 10- 


-3 


+ 1.98X 10- 


-4 


132.5 


10.00 


5.51 X 10- 


-4 


+ 8.06 X 10- 


-5 


95.5 


10.25 


2.83 X 10- 


-4 


+ 4.60 X 10- 


-5 


61.5 


10.50 


6.26 X 10 


-5 


± 1.59X 10- 


-5 


21.0 


10.75 


6.40 X 10 


-6 


± 4.66 X 10- 


-6 


3.1 


11.00 


8.17 X 10- 


-6 


± 5.23 X 10- 


-6 


3.0 



The entries are: Column 1: Center of luminosity bin. Column 2: 
Bivariate FIR luminosity function and its error. Column 3: Number 
of galaxies in the bin. The numbers are not integer due to the survival 
analysis applied. 

Rabaca ( I1994> earlier pointed out the difficulty with using 
a Schechter function to adequately describe nurture-affected 
samples. With only three galaxies above log(LFiR/L0) = 11.0 
our sample is well fit by a Schechter function. 

4.4. LpiR and Lb 

Figure0plots LpiR vs. Lg for the optically complete sample de- 
fined in Sect.EI] Lg was calculated as Lb = 10('-^5-o.4m„_„j 
in units of solar bolometric luminosity where Mzw-con is the 
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l = g(Lu) [L„] 



Fig. 7. LpiR vs. Lb for the optically complete, magnitude lim- 
ited subsample (n - 701, see Sect. 14.21 for exact definition). 
The full line indicates the best-fit bisector slope derived with 
ASURV, the dotted line shows the result of the regression 
adopting Lg, and the dashed line adopting LpiR (dashed) as in- 
dependent variable. 

absolute Zwicky magnitude coiTected for systematic errors. 
Galactic extinction, internal extinction and with K-coiTection 
applied (see Sect. 14. H and Paper I). This definition^ provides an 
estimate of the blue luminosity (vLy) at 4400 A. In Fig.Qwe 
see scatter due to measurement uncertainties and intrinsic dis- 
persion. The latter should be minimised as much as possible to 
nature-driven dispersion. In this sample we have reduced dis- 
persion due to both one-on-one interactions and environmental 
density. Most galaxies lie close to the coiTelation with a dis- 
persion of 0.23 for the detected galaxies. There are four clear 
outliers close to Lfir= 10"Lq. Three of them have been clas- 
sified as possibly interacting in Paper II. 

We applied linear regression analysis in order to estimate 
the functional relationship between the two variables. Since we 
are interested in the physical relation between the two vari- 
ables and since both variables have intrinsic uncertainties we 
decided to use a symmetric method. We derived the regression 
coefficients for both LpiR vs. Lb and Lb vs. Lfir using ASURV 
and calculated the bisector regression line shown in Fig.0from 
these, following the formula in Isobe et al. ( 1990l. We used the 
Schmitt's binning method (Schmitt ll985> as the only method 
offered by ASURV able to deal with censored data in the in- 
dependent variable. We note however, that for the cases where 
the other two methods (estimation-maximisation method, and 
Buckley-James method) could be applied, a satisfactory agree- 
ment was found. The results for the linear regression (ASURV 
bisector) are listed in Table |6] The alternate approach would 

- Note that this definition differs by a factor of 1 .7 from the defini- 
tion used in Paper I which was normalised to the solar luminosity in 
the blue. 
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Table 6. Correlation analysis of LpjR vs. Lg- 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Sample 


n 


"up 


slope 


intercept 


slope 


intercept 








(bisector) 


(bisector) 


(Lb indep.) 


{Lb indep.) 


Total 


701 


312 


1.41±0.02 


-4.55±0.25 


1.11+0.03 


-1.57+0.34 


Sa~Sab (1-2) 


49 


18 


1.37±0.09 


-4.29+0.87 


0.87+0.23 


0.76+1.67 


Sb-Sc (3-5) 


452 


181 


1.35±0.03 


-3.98+0.31 


1.04+0.06 


-0.77+0.57 


Scd-Im (6-10) 


115 


49 


1.25±0.03 


-2.99+0.30 


1.16+0.05 


-2.09+0.58 


Interacting 


14 


2 


1.52±0.12 


-5.25+1.15 


1.43+0.13 


-4.34+1.32 



The slope and intercept are defined as: log(LpiR) = log(Lg) x slope + intercept. The entries are: Column 1: Subsamples considered. All 
subsamples are selected from the optically complete, magnitude limited subsample (see Sect. l4.lt . In the early-type subsamples (E and SO) 
the relative number of undetected galaxies in LpjR is very high so that a regression slope could not be determined. Column 2: Total number of 
galaxies in the respective samples. Column 3: Number of galaxies with upper limits in FIR. Column 4: Bisector slope and its error of the best-fit 
regression line derived with the Schmitt binning method in the ASURV package. Column 5: Bisector intercept. Column 6: Slope and its error 
of the best-fit regression line derived with the Schmitt binning method in the ASURV package adopting Lb as independent variable. Column 7: 
Bisector intercept adopting Lg as independent variable. 



be to compute the regression assuming that optical luminos- 
ity is the independent variable. The results are also listed in 
Table |6] and show that the conclusions drawn in the follow- 
ing would not be substantially changed if Lb had been adopted 
as independent variable. The best-fit slope for the entire sam- 
ple is LpiR L^'*'^^""^. Our slope is shallower than the one 
found by Perea et al. til997j for a smaller subsample of the 
CIG, Lb oc Lpj^^^"^ (giving a slope of the inverse relation of 
LpiR ^B^^)- The main reason for this difference is our use of 
the bisector slope, whereas Perea et al. Cl997j derived the slope 
with LpiR as independent variable. With our larger sample we 
derive a similar slope when adopting LpiR as independent vari- 
able (Lfi oc Lpj^=^" "-'). For the present data set, however, we 
think that the bisector slope (or Lb as independent variable) is 
the better choice for investigating the functional relation be- 
tween both variables. A possible explanation for the slope >1, 
suggested by Perea et al. ( I1997> . is an increase of the dust ex- 
tinction with galaxy luminosity, yielding a faster increase of 
the FIR emission in comparison to the extinction-affected blue 
luminosity. An alternative reason could be the increase of the 
recent star formation (SF) activity (traced by Lfir) with galaxy 
luminosity. 

Figure |8] presents LpiR vs. Lb for subsamples of differ- 
ent Hubble types. Due to the low detection rate for early-type 
galaxies (E-SO), no reliable regression fit could be derived for 
this subsample. The correlation for the early-types shows ev- 
idence for a composite population with typical FIR deficient 
galaxies below the superimposed regression line and overlumi- 
nous galaxies, showing a roughly linear correlation spanning 
2 dex, above the line. As mentioned before, IR overluminous 
early-type galaxies must be regarded with suspicion until their 
morphologies are confirmed with higher resolution images than 
the POSS2 used for our morphology revision. At the same time, 
bona fide isolated early-types are of particular interest in view 
of ideas that see all of them as a product of nurture (merg- 
ers/stripping/harassment). There are only small differences in 
the measured slopes (see Table|6land Fig. |S} of least-squares 
regression lines as a function of Hubble type. 



0,5 - 




-1.5 - 



_2 I I I I I I I I I I I I I I I I I I 

-5 5 10 

Hubble type 

Fig. 9.R = log{Lfij^/LB), for the optically complete sample as a 
function of Hubble type. Only detections are shown. The open 
triangles give the mean value for each Hubble type, calculated 
with ASURV, taking the upper limits into account. The open 
squares are the median values for the detections. 

Finally, we derive the distribution and the mean value of 
R = log{LpiR/LB), a variable frequently used as an indicator of 
SF activity. In Table |4] we list the average and median values 
of R, together with those of LpiR and Lb derived for different 
subsamples. Figure|9lshows Rasa function of the morphologi- 
cal subtypes. No clear trend of < 7? > is found within the spiral 
galaxies with < R > essentially constant between T - 1-7 
(Sa-Sd). < R >, as well as < LpiR > (Fig.|5|l, is significantly 
lower for early-types (E and SO) although values derived using 
survival analysis might be uncertain due to the large number of 
upper limits. This means that early-type galaxies have a lack 
in FIR emission with respect to their blue luminosity with the 
ones showing values similar to spirals possibly being misclas- 
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Fig. 8. LpiR versus Lg for subsamples of different Hubble types. The fuU line is the bisector fit for the total AMIGA sample 
presented in Fig. whereas the dashed line is the bisector fit for the present subsample. For the early type subsample (E and 
E-SOa) no rehable regression fit could be derived due to the large number of upper limits. 



sified spirals. Late type galaxies (Sd-Im) are on average less 
luminous both in LpiR and Lb, but show the same < R > as 
spirals. 

4.5. IRAS colours 

IRAS flux ratios provide another potentially useful diagnostic. 
Feo/Fioo (Telesco et al..l988j. Fis/Feo (XS91) and F12/F25 
(Bushouse et al. I1988> have been used as environmental diag- 
nostics. For example, Feo/^ioo measures the dust temperature 
and has been found to increase with the level of star formation 
activity (de Jong et al. .l984j . /^is/^'eo is an indicator for AGN 



activity with values above ^25/^50 - 0.3 regarded as indicative 
of a Seyfert nucleus (de Grijp et al. ll985t . 

Figure ^1 presents histograms of different IRAS colours 
for our optically complete subsample. The average and median 
values are listed in Table The flux ratios log(F6o/^^ioo) and 
log(Fi2/f25) show a relatively symmetric distribution around 
the peak values. On the other hand, log(F25/f eo) exhibits a tail 
towards high values. The relative intensity of this tail weak- 
ens when only including detections with a higher S/N (we 
used S/N > 7 as a test), suggesting that part of it might 
be due to uncertain values, mainly at 25 fim. Another pos- 
sible reason for high values of ^25/^^60 can be the presence 
of AGNs, following the finding of de Grijp et al. J1985> that 
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Table 7. Mean and median values of IRAS colours. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Sample 






< log(F6o/-Fioo) > 
med(log(F6o/Fioo)) 




< log(F25/F6o) > 

med(log(F25/F6o)) 




"up 


< log(Fi2/F25) > 

med(log(Fi2/F25)) 


Total 


468 


76 


-0.42+0.01 
-0.45 


343 


-0.87+ 0.02 
-0.83 


126 


67 


-0.33+0.03 
-0.18 


S/Im (T = 1-10) 


443 


72 


-0.43±0.01 
-0.45 


326 


-0.87+0.02 
-0.83 


118 


65 


-0.33+0.03 
-0.16 


E(-5) 


9 


3 


-0.23+0.06 
-0.23 


4 


-0.73+0.09 
-0.74 


5 


2 


-0.47+0.07 
-0.48 


SO (-2) 


10 


1 


-0.39±0.06 
-0.35 


8 


-1.02+0.004 
-1.02 


2 





-0.16+0.01 
-0.16 


SOa (0) 


6 





-0.36±0.07 
-0.27 


5 


-0.98* 
-0.98 


1 





-0.27* 
-0.27 


Sa(l) 


9 


2 


-0.43±0.04 
-0.42 


7 


-0.79+0.05 
0.71 


2 


1 


-0.32* 
-0.32 


Sab (2) 


27 


3 


-0.42±0.03 
-0.45 


18 


-0.81+0.04 
-0.72 


9 


4 


-0.39+0.09 
-0.22 


Sb (3) 


88 


15 


-0.41+0.02 
-0.44 


63 


-0.87+0.03 
-0.89 


25 


17 


-0.50+0.10 
-0.18 


Sbc (4) 


104 


17 


-0.45+0.01 
-0.46 


77 


-0.83+0.03 
-0.78 


27 


15 


-0.27+0.04 
-0.15 


Sc (5) 


138 


24 


-0.46+0.01 
-0.48 


107 


-0.87+0.03 
-0.88 


32 


14 


-0.24+0.05 
-0.10 


Scd (6) 


34 


6 


-0.45+0.02 
-0.44 


24 


-0.87+0.06 
-0.82 


10 


6 


-0.36+0.08 
-0.17 


Sd (7) 


21 


2 


-0.39+0.03 
-0.40 


13 


-0.85+0.07 
-0.79 


8 


5 


-0.23+0.003 
-0.23 


Sdm (8) 


4 


1 


-0.40+0.02 
-0.40 


3 


-0.70+0.15 
-0.34 


1 


1 




Sm (9) 


5 


1 


-0.27+0.02 
-0.28 


3 


-1.19+0.08 
-1.24 


2 


1 


-0.28+0.04 
-0.22 


Im(lO) 


13 


1 


-0.31+0.03 
-0.35 


11 


-0.95+0.08 
-0.80 


2 


1 


-0.09+0.10 
0.06 


Interacting 


14 


2 


-0.36+0.03 
-0.39 


10 


-0.87+ 0.03 
-0.89 


4 





-0.32+0.08 
-0.34 



For the entries marked with "*" ASURV was not able to calculate an error. A "_" means that the entry could not be calculated due to the 
low number of detections. For ratios involving Fgo, only galaxies with detections at 60 fjm are taken into account and for log(Fi2/F25) only 
galaxies with detections at 25 fim. The entries are: Column 1: Considered subsample. All subsamples are selected from the optically complete, 
magnitude limited subsample. The interacting subsample consists of galaxies excluded from the CIG in Paper II (see Sect. 14. U . Column 2 and 
7: Total number of galaxies in the subsample. Column 3, 5, and 8: Number of galaxies with upper limits. Column 4, 6 and 9: First row: Mean 
value of the ratio, using the Kaplan-Maier estimator from ASURV. Second row: Median value of the same ratio, only for detections. 



galaxies with ^25/^60 > 0.3 are very likely to host an AGN. 
We have checked the values of F25/F60 for galaxies with an 
AGN listed in Sabater et al. (in prep.). Their list includes galax- 
ies catalogued to have an AGN in NED or in the Veron-Cetty 
Catalogue of Quasars and Active Nuclei (Veron-Cetty & Veron 
[2003 1, as well as radio-excess galaxies with radio luminosities 
more than 5 times the values predicted by the radio-FIR cor- 
relation and which are likely to be radio-loud quasars (Sopp 
& Alexander 1991 ). We found that 10 out of 11 active galax- 
ies with detections at both 25 fim and 60 fim have values of 
log(F25/F6o) > -0.7, the value where the departure from sym- 
metry in the distribution of F25/F(,o starts to be noticeable. 
Furthermore, 10 out of 14 galaxies with upper limits at F25 
might lie above this threshold, but the upper limit at F25 makes 
a firm conclusion impossible. Thus, even though the absolute 
number of galaxies with known AGNs is not enough to explain 



the tail towards high ^25/^60, they might be responsible for 
part of it. 

In Fig. ^2 we show the different IRAS colours as a function 
of Hubble type and in Table {7\ we list the average and mean 
values. We notice the following: 

- The range in log(F6o/^'ioo) occupied by most galaxies 
is quite narrow, with almost all objects having -0.7 ^ 
log(F6o/i^ioo) 5 -0.2. 

- For log{F(,Q/Fio()), we find a clear trend with Hubble type. 
The value is highest for ellipticals (type -5), decreasing 
towards spirals and increasing again for late-type galaxies, 
starting from type 7-8, until irregulars (type 10). 

- There is no significant trend in Hubble type visible for 
log(F25/^6o) or log(Fi2/F25). The low number of detec- 
tions might be the reason. We neither found a trend for 
\og(Fi2/F^o) (not shown here), for which we derived a 



Lisenfeld et al.: The AMIGA sample of isolated galaxies III 



13 



mean value for the total sample of -1.13+0.02, and very 
similar values for each Hubble type individually. 

5. Discussion 

5.1. Comparison to other non-interacting samples 
5.1 .1 . LpiR and Lb 

We compare the distribution of the FIR luminosity and of R to 
that of the galaxy sample of the Center of Astrophysics (CfA, 
Huchra et al. ll983t . whose FIR properties, based on data of the 
IRAS FSC, were studied in Thuan & Sauvage ( fT992t (here- 
after TS92) and Sauvage & Thuan (1992i (hereafter ST92). 
The CfA sample consists of 2445 galaxies representing a com- 
plete flux-limited sample (mz„ < 14.5) selected in Galactic co- 
ordinates. No selection with respect to environment was car- 
ried out. In order to properly compare the two data sets we 
applied the same magnitude cutoff' as in TS92 (otzw < 14.5, in 
uncorrected Zwicky magnitudes), to our sample. We then com- 
pared the velocity distribution of these two samples (the CfA 
sample and our adapted sample) and found a very good agree- 
ment, with only two differences: in the CfA sample the peak 
at ~5000 kms"' is missing due to their restriction in coordi- 
nates which avoids the region of the Perseus-Pisces superclus- 
ter responsible for this peak. Furthermore, in our sample with 
the above magnitude restriction there were no galaxies beyond 
8500 km s"', whereas about 4% of the galaxies in the CfA sam- 
ple have velocities above this value. We checked the effect of 
excluding these high velocity galaxies in the CfA sample on 
the subsequent results and did not notice any significant differ- 
ences. 

In order to correctly compare the luminosity distributions, 
we derived the distances for the CfA sample in the same way 
as for our galaxies: for close-by galaxies (Vhei < 1000 kms ') 
(for which we used redshift-independent distances from the lit- 
erature) we adopted the distances given by TS91, who used a 
Virgo-infall model to calculate them. For galaxies with Vhei > 
1000 kms"' we calculated the velocities after the 3K correc- 
tion, VsK, in the same way as for the AMIGA sample (see 
Paper I), and derived the distances as D = V^k/Hq- We used 
the same Hubble constant {Ho - 75 kms^'Mpc"') in both 
samples. 

As a test to find possible systematic differences we com- 
pared the distances, Lfir, Lb and R for those galaxies that are 
common in both samples (total: n - 98, with FIR detections in 
both samples: n - 87). TS92 used to derive Lb- For the 
CfA sample we estimated the corrected Zwicky magnitudes 
from using the linear relation found between both quanti- 
ties in paper I. Then we calculated Lb with the same formula 
as for the AMIGA sample. The calculation of LpiR was also 
done in the same way for both samples. For the 98 galaxies we 
found an excellent correlation between the distances used by 
us and those based on data of TS92, with a correlation coeffi- 
cient of 0.995 and a slope of 1.01 + 0.01. We also found a very 
good correlation between our values of Lfir and the values de- 
rived by TS92 (correlation coefficient of 0.96 for detections) as 
well as for Lb (correlation coefficient of 0.90) and for R (cor- 
relation coefficient of 0.85 for detections). The mean values of 



log(LFiR), logiLs) and R for the galaxies in common practically 
agree (see Table|8}, showing that a comparison of both data sets 
is justified. 

In Fig. El we show a comparison of our distribution of 
log(LFiR) to that of the CfA sample. Above log(LFiR/Lo) - 
10.2 a clear excess of CfA galaxies in comparison to our sam- 
ple is visible. In Table |H1 we list the mean values. The dif- 
ference between the mean value of log(LFiR) of the AMIGA 
and the CfA sample is 0.21-0.26 (with and without taking into 
account upper limits) which is a difference of 3^cr. We per- 
formed statistical two-sample tests in the package ASURV and 
found that the two distributions were different with a proba- 
bility between 97.22% (Logrank test) and 99.87% (Gehan's 
GeneraUsed Wilcoxon Test). The maximum probabiUty in- 
creases to > 99.95 % when only detections are taken into ac- 
count. We also performed a Kolmogoroff-Smirnoff test on the 
detected data points and derived a probability of more than 
99.75% that the mean values of Lfir are different. Therefore, 
there is strong statistical evidence that the AMIGA sample has 
a lower Lfir than the CfA sample, which is comparable but not 
selected with respect to the environment. This suggests that the 
FIR luminosity is a variable driven by interaction. 

The comparison of the distribution of R is shown in Fig. 1131 
We notice that the mean value of R is higher for the CfA sample 
than for the AMIGA sample. The difference is 0. 12-0. 14 (with 
and without upper limits) (see Table |8j which corresponds to 
4-7 cr. This difference has its origin in the higher value for Lfir 
of the CfA sample, as the mean values for Lb are very simi- 
lar (see Table |8j, and the distribution of Lb for both samples 
(not shown here) practically agrees. We performed the statisti- 
cal two-sample tests in the package ASURV and found that the 
two distributions of R were different with a probability between 
97.3% (Logrank test) and > 99.95% (Gehan's Generalised 
Wilcoxon Test). On the other hand, the tests showed that the 
distributions of Lb were identical with a non-negligible proba- 
bility (28-68%) confirming that the difference in R has its ori- 
gin in Lfir. Performing a Kolmogoroff-Smirnoff test on the de- 
tected data confirmed these results, yielding a probability of 
more than 99.999% that the mean values of R are different. 

5.1.2. IRAS colours 

We compared the distribution of the IRAS colours to the re- 
sults found for the IRAS Bright Galaxy Sample (BGS, Sanders 
et al. 2003 1. The value oflog(Fa)/Fmo) of the BGS is higher by 
about 0.2-0.3 with respect to our sample. This is not surprising, 
since the BGS contains galaxies in a more active star forming 
phase than the CIGs. The peaks of the distribution of the other 
colours in the BGS are very similar to ours, the only excep- 
tion being the asymmetric tail towards high values of F 25 1 Fa) 
which is absent in the BGS. 

A comparison to the results of XS91 for a smaller subsam- 
ple of CIG galaxies (see Table|9j shows an excellent agreement 
of the values for \og{F I F ioq) and log(Fi2/F25)- Our value for 
log(F25/L'6o) is however slightly higher than that of XS91. We 
also compared our results to the CfA sample studied by ST92. 
For this aim, we produced again a different subsample, carry- 
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Table 8. Comparison to the CfA sample (Thuan & Sauvage [l992> . 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Subsample 


n 


< log(LB) > 




"up 


< log(LFIR) > 


O" FIR 


< \og(R) > 


o-R 


Total subsamples 


















AMIGA (all) 


207 


9.80±0.05 




28 


9.16±0.09 




-0.56+0.03 




AMIGA (del.) 


179 


9.83±0.04 


0.56 





9.38+0.05 


0.71 


-0.45+0.02 


0.31 


CfA (TS92) (all) 


1544 


9.89±0.01 




210 


9.42+0.02 




-0.44+0.02 




CfA (TS92) (det.) 


1334 


9.89±0.01 


0.54 





9.59+0.02 


0.73 


-0.31+0.01 


0.41 


Galaxies in common 


















AMIGA(all) 


98 


9.82±0.06 




4 


9.37+0.07 




-0.43+0.03 




AMIGA (det.) 


87 


9.84±0.06 


0.55 





9.44+0.07 


0.66 


-0.40+0.03 


0.29 


CfA (TS92) (all) 


98 


9.85±0.05 




7 


9.37+0.07 




-0.47+0.04 




CfA (TS92) (det.) 


87 


9.85±0.05 


0.50 





9.43+0.07 


0.62 


-0.42+0.04 


0.33 



The entries are: Column T. Sample considered. Both samples are selected with the same magnitude limit of (uncorrected) ra^^ < 14.5. The 
distances of the galaxies of the CfA sample have been derived in the same way as for the AMIGA sample (see Sect. 15.1.11 . We give both the 
results obtained with ASURV (first row) and the results with detections only (second row). Column 2: Total number of galaxies. Column 3: 
Mean value of log(LB) and its error. Column 4: Standard deviation of log(LB). Column 5: Number of galaxies with upper limits in Lfir. Column 
6: Mean value of log(LFiR) and its error. Column 7: Standard deviation of log(LFiR). Column 8: Mean value of R = log(LFiR/LB) and its error. 
Column 9: Standard deviation of log(LFiR/LB). 



Table 9. Comparison of IRAS colours to other studies. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Sample 


"/"up 


< log(^) > 


n/n„p 




"/"up 


<log(f^)> 


Isolated samples 














AMIGA total 


468/76 


-0.42+0.01 


468/343 


-0.87+0.02 


126/67 


-0.33+0.03 


AMIGA, only det. 


392/0 


-0.44+0.01 


125/0 


-0.76+0.02 


59/0 


-0.18+0.03 


XS91 CIG 


261/. 


-0.42+0.01 




-0.96+0.02 




-0.32+0.04 


Bushouse et al. )1988> isolated 


68/0 


-0.39 






34/0 


-0.21 


AMIGA m™(„„corr) < 14.5 


183/4 


-0.43+0.01 


183/90 


-0.88+0.02 


93/35 


-0.26+0.03 


AMIGA m^„(„„eo,r) < 14.5, (det.) 


179/0 


-0.43+0.01 


93/0 


-0.82+0.02 


58/0 


-0.18+0.02 


CfA (ST92) 


1465/131 


-0.42+0.004 


1465/771 


-0.94+0.01 


706/154 


-0.26+0.01 


CfA (ST92), (det) 


1330/0 


-0.43+0.004 


694/0 


-0.87+0.01 


552/0 


-0.21+0.01 


Interacting samples 














AMIGA interacting 


14/2 


-0.36+0.03 


14/10 


-0.87+0.03 


4/0 


-0.32+0.08 


XS91 wide pairs 




-0.39+0.01 




-0.93+0.02 




-0.52+0.07 


XS91 close pairs 




-0.31+0.01 




-0.93+0.02 




-0.65+0.08 


Bushouse et al. ( 19881 inter. 


98/0 


-0.33 






48/0 


-0.42 



The entries are: Column I: Sample considered. Columns 2, 4, 6: Total number of galaxies and number of galaxies with upper limits. Columns 3, 
5, 7: IRAS colour. For ratios involving Fso, only data points with detections at this wavelength were considered, whereas in < log(Fi2/F25) > 
only data points with detections at 25 yum were taken into account. A "_" means that the corresponding data point was not given in the reference. 



ing out the same magnitude cut (< 14.5 in uncorrected Zwicky 
magnitude) as in ST92. We found a very good agreement for 
log(F6o/f lOo) (see Table|9jl, and for log(Fi2//^25)- With respect 
to log(F25/F6o), we derived a slightly higher value for our sam- 
ple. However, we consider the significance of this difference 
low due to the large number of upper limits. 

ST92 found in their analysis the same trend with Hubble 
type for log(F6o/f too) as we did. The value that they found 
for ellipticals, log{Fa)/Fioo) = -0.38, is slightly higher than 
ours whereas their value for irregulars, log(F6o//^ioo) = -0.32, 
agrees very well. Also their values for spirals (between -0.45 
and -0.47 for T = 2-5) are very close to ours. They explained 
the high Feo/^ioo ratio in ellipticals by the concentration of the 
dust in the central regions where the radiation field is high, pro- 
ducing in this way a higher dust temperature. A high Feo/^^ioo 
ratio for irregulars has been found in other studies as well (e.g. 



Melisse & Israel i 1994t and can be understood as a lack of "cir- 
rus" emission with respect to FIR emission from H ii regions. 

5.2. Comparison to interacting galaxies 
5.2.1 . LpiR and Lg 

One of the motivations for refining and studying a large sample 
of isolated galaxies is to better define a baseline against which 
effects of environment could be quantified. Both mean IR di- 
agnostic measures and their dispersion are of interest in this 
context. AMIGA began with a CIG sample selected to avoid 
as much as possible near neighbours. Yet visual reevaluation 
of the optical morphologies for the sample using P0SS2/SDSS 
revealed 32 objects with clear signs of interaction (Paper II). 
These galaxies have been excluded from the present effort to 



Lisenfeld et al.: The AMIGA sample of isolated galaxies III 



15 



characterise the isolated sample but offer a useful internal com- 
parison sample to measure the sensitivity of the IR diagnos- 
tics to environment. Fourteen of these 32 galaxies have IRAS 
measures and apparent magnitudes between 1 1 and 15. Table|3 
shows that the mean logCL^) is almost identical to the isolated 
sample while log(LFiR) is brighter by ~0.7 (respectively 0.6 
when comparing to the spiral/irregular subsample). Similarly 
the mean FlR-to-optical flux ratio, < R >,is 0.50 (respectively 
0.43 for the spiral/irregular subsample) higher This shows that 
there is significant difference in R between isolated and inter- 
acting galaxies due to an enhancement in Lfir of the latter. 

Figure shows the correlation between LpiR and Lg for 
this interacting subsample. Regression analysis yields a steeper 
slope (see Table |5Jl, as is usually found for interacting sam- 
ples, indicating that LpiR increases faster as a function of Lg in 
comparison to the non-interacting sample. The reason for this 
FIR excess is most likely an enhancement of (dust-enshrouded) 
SF in interacting galaxies. This is consistent with the results in 
Perea et al. ( 1 997) where a FIR enhancement was found for 
a sample of perturbed galaxies. We included their weakly and 
strongly interacting samples in Fig.O The effect is strongest 
for their strongly perturbed sample. The average FIR excess 
(i.e. average deviation from the AMIGA regression line) for 
the strong and weak interacting samples are 0.87 and 0.49 re- 
spectively. 

XS91 compared mean FIR-to-optical flux ratios and found 
a much smaller, but significant, difference between spiral-spiral 
pairs {R - -0.17) and a late-type subsample from the CIG 
(R = -0.30). The difference increased when considering only 
close pairs (R - -0.02). In order to compare their result 
to ours, we have to take into account that they used uncor- 
rected Zwicky magnitudes. The correction that we applied (see 
Paper I) changed the Zwicky magnitudes by on average -0.67 
magnitudes, corresponding to a change in R of -0.27 dex. 
Taking this into account, the value R for the CIG subsample 
of XS91 is practically the same ours, whereas R for the pair 
sample in XS91 is below (~0. 2-0.4) the value of our interact- 
ing sample. Due to the small size of our interacting sample and 
the different selection (the study of XS91 restricted the envi- 
ronmental signature to the effects of one-on-one interactions) 
we do not want to draw any conclusions from this difference. 

5.2.2. IRAS colours 

We found a slightly higher value of log(F6o//^ioo) for the 
possibly interacting galaxies in the CIG than for the total 
AMIGA sample (see Table0. The difference is, however, only 
2o", and thus not statistically significant. In the other colours 
(log(F25/f 6()), log(f i2//^25) or logiFn/Fa))) we found, within 
the errors, no difference between both samples. 

A trend of higher log(F6o//^ioo) values in interacting galax- 
ies has been found in other studies (see Table |9jl. XS91 com- 
pared thek subsample of CIGs to paired galaxies. They found 
a value very close to our interacting sample for wide late-type 
pairs (i.e. distance between partners larger than the diameter 
of the primary) and an even higher value, significantly higher 
than for the value for the CIGs, for close late-type pairs (i.e. 



Table 10. Galaxies with log(LFiR/io) > 10.5 



(1) 


(2) 


(3) 




(4) 


CIG 


log(LFiR/Lo) 


Hubble type 


Comment 


55 


11.12 


Sc (T = 


5) 


I/A + Sy/LINER 


143 


10.86 


She (T = 


:4) 


I/A? lopsided spiral 


148 


10.56 


Sbc (T = 


:4) 


I/A 


232 


10.69 


Sc {T = 


5) 


I/A? 


302 


11.02 


Sc (T = 


5) 


I/A? peculiar 


361 


10.51 


Sab {T = 


■■2) 


Isolated spiral? 


444 


10.54 


Sb (T = 


3) 


Possibly Sm companion 


709 


10.53 


Sc (T = 


5) 


Sm companion nearby 


715 


10.76 


Sc (T = 


5) 


I/A 


829 


10.51 


Sb{T = 


3) 


Blue Compact 


841 


10.58 


SO (r = 


-2) 


large inclined SO 


866 


10.95 


Sb (T = 


3) 


Isolated, LINER 



The entries are: Column 1: CIG number. Column 2: FIR luminosity. 
Column 3: Hubble type. Column 4: Comment after visual inspection 
of optical images, and consultation of NED. I/A means interacting, 
I/A?: possibly interacting, Sy: Seyfert galaxy 



of distance between partners less than the diameter of the pri- 
mary). Bushouse et al. ( 1988) studied the MIR/FIR properties 
of a sample of 109 colliding galaxies and compared it to more 
isolated galaxies from the sample of Kennicutt & Kent J1983> . 
In their study, they only took into account IRAS detections. 
Their values for log(F6o/Fioo), both for the interacting as well 
as for the comparison sample, is in reasonable agreement with 
our study and that of XS91. 

We did not find a significant difference in neither ^25/^^60 
nor F12/F25 between isolated and interacting galaxies. In con- 
trast to this, XS91 obtained a lower value for F12/F25 both 
for the close and the wide pair samples. Also Bushouse et al. 
( 1988 1 found a higher value for F12/F25 in the interacting sam- 
ple, although their result has to be taken with caution because 
only detections were included (thereby skewing the result to 
higher values). Given the very low number of galaxies with de- 
tections in our interacting subsample, the significance of these 
differences is, however, not statistically meaningful. 

5.3. Nature of the FIR brightest galaxies 

MIR/FIR measures have been found to be sensitive diagnos- 
tics of enhanced SF. Since SF can be greatly enhanced by the 
presence of companions we can ask if any of the most FIR lu- 
minous galaxies in our sample are really isolated. Naturally, as 
discussed in Paper II, we are limited in our ability to describe 
galaxy form and to detect close companions by the quality of 
the available images. Table [Tol lists the FIR brightest galaxies 
(with log(LFiR/i'o) > 10.5) in our sample and the results of re- 
examination of POSS2 and other data sources for them. More 
than half of this sample are likely to involve interacting sys- 
tems. A very few may represent the most luminous examples of 
self-induced star formation or the IR measures are dominated 
by a yet undetected active nucleus (many do not have published 
spectra). CIG 709 epitomises another issue raised in Paper II. 
While it appears to be isolated from similarly sized objects, it 
shows striking structural asymmetries and a small companion 
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about one diameter distant. It is yet unclear whether such small 
companions are capable of significantly enhancing star forma- 
tion and producing structural deformations in massive spirals 
(see also Espada et al. 2005 1. 

We furthermore inspected the most extreme outliers from 
the regression fit between Lfir and Lb- The outliers on the high 
side come under immediate suspicion as interactors that were 
missed in the morphology survey. Only two show log(LFiR) > 
1 1 .0 and both are possible interacting systems - in NED they 
are described as Seyfert/LINER (CIG 55) and Hii galaxy 
(CIG 302). Several others are candidate interactors or have an 
active nucleus. Galaxies falling below the regression are either 
early-types (E and SO), known to be deficient in FIR emission 
or highly inclined galaxies. The internal extinction correction 
for such objects is large and uncertain suggesting that the low 
Lpir/Lb measures may be due to overcorrection of the blue 
magnitude. 

6. Conclusions 

We present a MlR/FlR analysis of a sample of the most isolated 
galaxies in the local universe, obtained from the Catalogue of 
Isolated Galaxies (CIG). This analysis is part of our multi- 
wavelength study of the properties of the interstellar medium 
of this sample and involves ADDSCAN/SCANPI reprocessing 
of IRAS data for all 1030 galaxies (out of 1050 in the CIG) 
covered by IRAS. We increased the detection rate with respect 
to the PSC and FSC in all IRAS bands and present our AMIGA 
sample of 701 CIG galaxies as the best available control sam- 
ple against which to evaluate the IR signatures of environment 
in local galaxy samples. Our sample is large enough to permit 
comparison of IR properties for galaxy morphology subclasses. 
Our main results are the following: 

1. The galaxies in our sample have modest FIR luminosi- 
ties, with only 14 objects (corresponding to <2% of the 
sample) above log(LFiR/i'0) = 10.5. The mean log(LFiR) 
of our sample is 3-4cr (0.21-0.26 dex) below the corre- 
sponding value for the CfA sample studied by TS92 and 
ST92, which is a sample of nearby galaxies similar to ours, 
but selected without considering the environment. In addi- 
tion, a lower value (by 0.12-0.14 dex, corresponding to 4- 
7cr), compared to the CfA sample, was found for the mean 
R = log(Lpif!_/LB) of our sample. This suggests that the FIR 
emission is a parameter driven by interaction and that our 
sample of isolated galaxies shows a value close to the low- 
est possible. 

2. We find evidence for a systematic increase in FIR lumi- 
nosity from type SO/a to Sc followed by a decline for later 
types (dominated by dwarf galaxies), possibly reflecting 
lower dust masses in those galaxies or less eflicient star 
formation. At the same time, R is essentially constant for 
all Hubble types later than SO/a suggesting that the SF efli- 
ciency in isolated spirals and irregulars is roughly constant. 

3. Early-type galaxies (E and SO) show a lower average R than 
the spirals. We can divide them into two populations: 1) un- 
detected galaxies (the majority) with low upper limits in R 
and 2) a population of early type with Lfir and R that are 



similar to spiral galaxies. The latter galaxies require con- 
firmation of the assigned early-types. If real, they represent 
an interesting class of isolated galaxies. 

4. We calculated the bivariate FIR luminosity function which 
was found to be in good agreement with previous stud- 
ies (XS91) based on a smaller subsample of the CIG. The 
FIRLF is dominated by moderately FIR luminous galaxies 
(only 3 objects have log(LFiR/L0) > 11.0) and is well de- 
scribed by a Schechter function. This contrasts to results 
obtained for FIR selected samples (e.g. Soifer et al. 1986, 
Sanders et al. l2003> where a double power law is needed to 
achieve a fit to the high-luminosity end of the FIRLF. 

5. We found a correlation between Lfir and Lb with a slope 
above 1 (Lfir °^ Lb^ '^^) with only modest variations as a 
function of Hubble type. Possible reasons for the slope be- 
ing >1 are an increase in extinction or an enhancement of 
the recent SF activity with galaxy luminosity. 

6. We identified a small population of possibly interacting 
galaxies in the CIG (Paper II) and these show a significantly 
higher mean FIR luminosity than the rest of the sample. 
They lie above the regression line derived for our optically 
selected CIG sample. The same was found for samples of 
interacting galaxies from Perea et al. J1997> . This suggests 
that the FIR emission is enhanced due to the interacting. 

7. We found a trend of Feo/^^ioo with Hubble type: ellipti- 
cal galaxies and irregular galaxies have a higher /^eo/Z^ioo 
than spirals, indicating a hotter dust temperature. For the el- 
lipticals the most likely reason is the higher concentration 
of dust towards the inner regions of the galaxies (ST92) 
whereas in irregulars, a lack of cirrus emission is the most 
probable cause. 

8. The value of Feo/Fioo of the AMIGA sample was found to 
be lower than that of interacting samples from the literature 
(XS91 and Bushouse 1988 1 indicating that interaction can 
increase the dust temperature. 

As the largest and most isolated sample in the local 
Universe, AMIGA can serve as a valuable control when as- 
sessing the effects of environment on other local samples of 
galaxies. This can in turn clarify our interpretation of the FIR 
signature in samples at higher redshift. 
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Fig. 10. IRAS Colours for the AMIGA sample. The full line 

shows the histogram calculated with ASURV taking into ac- 
count censored data points, the grey-shaded area shows the 
detections and the dashed line the upper limits. Only galaxies 
with detections at 60 /im are considered for log(F6o/^^ioo) and 
log(F25/F6o) , and only galaxies with detections at 25 jum for 

log(Fi2/F25). 
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Fig. 12. The percentage FIR luminosity distribution for the FIR 
detections in the AMIGA sample restricted to (uncorrected) 
Wz„ < 14.5 (shaded area) and the corresponding distribution 
for the CfA sample (dotted line). 
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Fig. 11. IRAS colours as a function of Hubble type for the op- 
tically complete sample. For log(F6o/^ioo) and log(F25/^^6o) 
only galaxies with detections at 60 fim are taken into account, 
and for log(Fi2/f25) only galaxies with detections at 25 jum. 
Only detected galaxies are plotted (crosses). The open squares 
indicate the mean values from Table0 calulated with ASURV 
and taking into account censored data points. The open trian- 
gles indicate the median values based on detections only. When 
no error bar is given, this could not be calculated due to the low 
number of detections. 



Fig. 13. The percentage distribution of R - logiLpi^/LB) for 
the FIR detections in the AMIGA sample restricted to (uncor- 
rected) m^w < 14.5 (shaded area) and the corresponding distri- 
bution for the CfA sample (dotted line). 
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Fig. 14. LpiR vs. Lb for various interacting subsamples. The 
filled squares and arrows denote the galaxies from the CIG 
showing signs of interaction (see Paper II). The dashed line 
is the regression fit to this subsample. The triangles indicate 
strongly interacting galaxies, and the crosses weakly interact- 
ing galaxies from Perea et al. (I1997> . The full line is the fit to 
the total AMIGA sample (n = 701) of Fig.El 



